Renormalized Surface Charge Density for a Strongly Charged Plate in Asymmetric 
Electrolytes: Asymptotic Exact Results in Poisson Boltzmann Theory 
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The Poisson-Boltzmann equation for a strongly charged plate inside a generic charge-asymmetric 
electrolyte is solved using the method of asymptotic matching. Both near field and far field asymp- 
totic behaviors of the potential are systematically analyzed. Using these expansions, the renormal- 
ized surface charge density is obtained as an asymptotic series in terms of the bare surface charge 
density. 
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I. INTRODUCTION 

A charged object inside an electrolyte attracts ions 
of opposite charge and repels ions of like charge. The 
total electrostatic potential, due to both the external 
charges and the electrolyte, is described by the Poisson- 
Boltzmann (PB) equation [l[ at the mean field level. 
Study of PB is a significant problem in both physical 
and biological sciences. Even though it has been clear in 
recent years that the Poisson-Boltzmann theory is inad- 
equate for strongly charged systems 0-Hl , it remains a 
benchmark theory in interpreting experimental data on 
charged macromolecules, particles, and surfaces. Despite 
of its importance, our analytic understanding of PB is 
very limited, largely due to its nonlinear nature. Even 
for the simplest problem of one charged plate, exact so- 
lutions are known only for cases of 1 : —1 and 2 : — 1 
electrolytes [f|. 

Consider an electrolyte with concentration c+ of m va- 
lence positive ions and concentration c_ of — n valence 
negative ions. It is asymmetric ii m =/= n. Inside this 
electrolyte, and at mean field level, the average electro- 
static potential near a charged plate satisfies the cele- 
brated Poisson-Boltzmann equation (PB) fl7| : 



&j>"(z) +c-nqe n l 3q 't' {z) 



c + mq e 



-mPq<p(z) _ (\X) 



where z is the coordinate perpendicular to the plate, 
P = l/k B T, and q = 1.6 X 10~ 19 C the charge of an 
electron. Because the dielectric constant of the plate 
is much smaller than that of the solvent, the standard 
electrostatic boundary condition reduces to that of Neu- 
mann: 



ecp'(z ) 



(1.2) 
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where zq is the location of the plate, and e dielectric 
constant of the solvent, and a the surface charge density 
on the plate. It will become clear below why we choose 
to locate the plate not at the origin of coordinate system. 
Note that the overall charge neutrality of the electrolyte 
imposes the constraint mc + = nc— . 

Far away from the charged plate, the potential must 
decay to its bulk value, which is chosen to be zero. There- 
fore in the far field the PB can be linearized: 



A<f> + n 2 <f> = 0, (1.3) 



where 



f-2 
-DB 



-(c + m 2 + c-n 2 ), (1.4) 



with Idb the Debye length. It is the characteristic length 
scale over which the electrostatic potential decays in the 
bulk. The far field potential is therefore: 

0fl 



-kSz 



z — >■ oo, 



(1.5) 



with Sz = z — zq the distance from the plate to the field 
point, an is called the renormalized surface charge den- 
sity. It acquires this name because Eq. (| 1 . 5|) is the solu- 
tion to linear PB Eq. (|1.3[) with renormalized boundary 
condition: 



£0'(^o) = o-R- 



(1.6) 



Physically, is the effective surface charge density we 
infer from the far field asymptotics if we use the linear 
PB equation. Evidently, the value of an can only be de- 
termined by comparing the linear solution Eq. (|1.5[) with 
the solution to the full nonlinear PB Eq. (jl.lj) . Calcula- 
tion of the renormalized charge density an as a function 
of the bare surface charge density a is the problem that 
we shall tackle in this work. 

It is convenient to rescale coordinate by the Debye 
length Iob an d to rewrite PB in terms of the dimen- 
sionless potential "J = f3q4>: 
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*"(z) + 



= 0. 



(1.7) 
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This equation can be integrated once, with the constant 
of integration determined by the boundary condition at 
infinity: *(oo) = *'(oo) = 0: 



1 



m + n 



1 1 

n m 



0. (1. 



The Neumann boundary condition Eq. (|1.2[) reduces to 

(1.9) 



* {zq) = = 77, 



where 77 is the dimensionless bare surface charge density. 
The renormalizcd boundary condition Eq. (jl.6l) then re- 
duces to 



where 



*'(*o) = 77*, 



(1.10) 



(1.11) 



is the dimensionless version of the renormalizcd charge 
density. One of the main purpose of this work is to cal- 
culate the functions 77^(77) for general cases of asymmet- 
ric electrolytes. As we shall find, these functions only 
depends on the type of electrolyte (m : —71), but are in- 
dependent of other parameters such as temperature and 
dielectric constant of the solvent. Using Eq. (jl.9|) and 
Eq. p. lip , we can represent the renormalizcd charge den- 
sity ur in terms of the function 77^(77) as 



or 



7T VRiv) 

pq 



en ( /3qa 

pq \ e« 



(1.12) 



which is sufficient to determine the leading order behav- 
ior of far field asymptotics Eq. (jl.5p . 

For (777,77,) = (1,1), and (2,1) the exact solutions to 
PB Eq. (|l.lj) arc known. From these, it is seen ijr ap- 
proaches a finite constant 7;^(±oo) as the bare charge 
density approaches ±00. 

7? fl (±oo) = ±4 1 : -1, 

t]r(oo) = 6 2 : — 1 positive plate, (1.13) 

t?h(-oo) = -^75 2 : -1 negative plate. 

Therefore in PB theory, the potential in the bulk remains 
finite even when the bare charge density goes to infinity. 
The full functions 77^(77) are also known. For details, sec 
Table HVl also see reference @. 

That 77#(±oo) saturates is a general property of non- 
linear PB theory, independent of the shape of the charged 
objects. Already exhibited by the exact solution 

to one dimensional PBE within symmetric electrolyte 
which has been known for a long time, the significance 
of this property has only been appreciated recently 0- 
[Tlj . As one important consequence, it implies that the 
interaction between two charged colloids saturates when 



their charges scale up, with their separation fixed. This 
was seen in numerical solution of PB equation (l2T | of a 
charged sphere, where the terminology charge renormal- 
ization was first introduced. 

In this work, we shall solve PB for one plate problem 
in a generic m : —n asymmetric electrolyte. We shall 
first develop asymptotic expansions of the electrostatic 
potential both in near field, for a plate with infinite bare 
surface charge density, located at the origin of axis z = 
(Sec. [II]). From these expansions, it becomes clear that 
the potential does remain finite even if the bare surface 
charge density is infinity. We shall then develop asymp- 
totic expansions in the far field, of which Eq. (|1.5p is the 
leading order (Sec. IIIip . In Sec.HV| we shall match the far 
field expansion with the near field expansion, and deter- 
mines the renormalizcd surface charge density 77^ for the 
infinitely charged plate. Finally in Sec. [Vj we exploit- 
ing the properties of near field expansion to obtain an 
asymptotic expansion of the renormalized charge density 
r/u in terms of the bare density 77. It is important to point 
out that our method can be immediately generalized to 
the cases of electrolyte mixtures. 

Before starting to solve Eq. (|1.8p . let us discuss two 
symmetries of this equation. Let us use <]7^/ - " to denote 
a solution to Eq. (|1.8p inside m : — n electrolyte, gener- 
ated by a positively /negatively charged plate. Eq. (|1.8[) 
is a first order ODE, and therefore contains only one pa- 
rameter in its general solution, which can be chosen to 
be the location of the plate. Therefore the solution to 
Eq. (|1.8p is unique up to the translation of the plate. 
Firstly, it is easy to check that — \]>™ ,_ " is a solution to 
PB inside an 77 : —m electrolyte, generated by a neg- 
atively/positively charged plate. We therefore find the 
following relation: 

^- m (z) = -*+'-™(z). (1.14) 

From this we obtain the following relation: 

V%'- n (v) = -VR- m (~V)- (1-15) 



As a consequence, we only need to solve Eq. (jl.8p for all 
cases with m > n. 

Secondly, if m and n have a common integer factor p, 
such that m = p 771, n = p n, then we can define \t = p^ ', 
which satisfies 

(1.16) 



777 + 77 



e"*-e-™* =0. 



Hence 



\z)=p-^f- n {z). 



(1.17) 



From this, we obtain the following relation between 
renormalized charge densities for the renormalized sur- 
face charge density: 

rf£- n {ri) - p-'vf-Hpvl (1-18) 
VR- n (±oo) = p- 1 V f- h (± o). (1.19) 

As a consequence, we only need to solve Eq. (|1.8p for all 
777, 77 that are relatively prime. 
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II. NEAR FIELD ASYMPTOTICS 

Sufficiently close to a strongly charged plate, the co- 
ions which carry like charges as the plate are strongly 
repelled so that their density is negligibly small. Con- 
sequently, we may ignore the corresponding term in the 
Poisson-Boltzmann equation. For a positively charged 
plate, the equation Eq. (|1.7p (with positively charged ions 
ignored) reduces to 



*o+(*)- 



1 



0, 



(2.1) 



which has a solution: 
*o+(*) = - 



2, 1, 2(m + n) , . 

— ]nz + -ln-i -. (2.2) 

n n n 



For a negatively charged plate, Eq. (|1.7[) reduces to 



1 



m + n 



0. 



(2.3) 



presence of counter-ions only. However, since the co-ions 
arc strongly repelled by the plate anyway, their existence 
is not going to change the nature of the singularity at 
3 = 0. The fact that the potential generated by an in- 
finitely positive/negative charged surface inside an elec- 
trolyte is finite for any z > is probably the most im- 
portant property of Poisson-Boltzmann theory. It was 
named (somewhat vag uely) as charge renormalization by 
Alexander et. al. [12J . This property of course can only 
be approximately applicable in real systems. All ions 
have finite volumes and the density of screening ions must 
saturate at the corresponding close packing value. For 
example, inside a ImM monovalent salt and near a sur- 
face with potential 0.3U, the ionic concentration would 
be about 10 5 mM according to the Poisson Boltzmann 
theory, which already exceeds the close packing limit, if 
every ion has (hydrated) volume of 30 A 3 . The Poisson- 
Boltzmann theory can be modified to incorporate finite 
size effects of ions [l3|, EH ■ 



which has a solution: 



. , 2 1 2(m + n) , . 

#o-(z) = — In 2 In— — ^— L . (2.4) 

m m m 

The most salient feature of these solutions are their 
singularities at z = 0. The potential diverges logarith- 
mically, while ion density diverges as z~ 2 as z — > 0. 
Therefore Eqs. (|2.2I2.4[) are the potentials generated by 
an infinitely charged plate located at the origin in the 



From now on we shall fix the infinitely charged plate 
at the origin. Eqs. (|2.2I2.4[) are the leading terms of sys- 
tematic near field expansions for the solutions to the full 
PB equation. Defining 



*±(z) = *o±(*) + M*) 



(2.5) 



and substitute these back into the full PB equation 
Eq. (|1.8[) . we obtain the equations for 9±(z): 



J 



2 mz 



1 



m + n 
1 



m + n 



IJrn + n) e+ _ _ 1 / 
n z z z m \ 



1 



n \ 2(m + n) 



r 



2(m + n) 
g _ 2(m + n) 



z " e 



1 1 



{e- ma - - 1) - - - — 
n m 



= 0, (2.6a) 
= 0. (2.6b) 



Using method of dominant balance |l5|, it is easy to see 
that for small z the leading terms for 9± scale as z 2 . How- 
ever, it would be incorrect to conclude that 9± can be 
expanded into Taylor series of z. For non-integer 2m jn 
(2n/m), the equation for 6 + (0-) contains nonanalytic 
term z 2m / n (z 2 ™/ m ), and hence the corresponding expan- 
sion of 9+ (9- ) must also contain this non-analytic term, 
as well as its integer powers. The explicit forms of the 
near field expansion for various cases of electrolytes are 
displayed in Table HI Using Wolfram Mathematica and 
substituting these expansions into the equations for 9±, 
Eqs. (|2.6[) , all coefficients are determined in terms of 02, 



which in turn can be determined exactly: 



1 

6m ' 



1 



(2.7) 



The near field expansion therefore does not contain any 
undetermined parameter. We summarize the coefficients 
in the Appendix [Cj 



III. FAR FIELD ASYMPTOTICS 

The simple form Eq. (|1.5[) is only the leading order 
term in the far field expansion. To find all higher order 
terms, it is convenient to introduce the following vari- 
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Electrolytes 


Plate 


Near field expansion of 0± 




3 : -1 


+ 


e+ = E*=i a ^ 2k 

Q- =E fe=0 «2+2fc/3^ +2fc/3 




4 : -1 


+ 


y-i x — > 94-k/ 1 ? 
8- =E fe= a 2+ fe /22 2+fe/2 




3 : -2 


+ 


0+ - Z^fe=2 fl fc 2 + Z Z^fe=l °2fcZ 








/i 2k , 4/3 i 2k 
°~ = Efc=l a 2fcZ + «' Efc=l»2fcZ 


8/3 2fc 


4 : -3 


+ 


= E fc =i oak* 3 * + * 8/3 E fc =i k^* 


+ ^ 16/3 E fc . 1 C2^ 2fe 






= Efc=i <>2kz 2k + z 3/2 E fe= i b 2k z 2h 


+ * 3 E fe=1 C2^ 2,: + 2 9/2 E fc = 1 d2^ 2fe 



TABLE I: Near field expansion of 0± as defined in Eq. (|2.5|) . where ^±(2) are defined in Eqs. (|2.2I2.4[) . Note the terms with 
fractional powers of z in the expansions. All coefficients are determined by the lowest order one aa, which is shown in Eq. (|2.7I) . 



ables: 

s = e~\ T(s) = e* (z) . (3.1) 
Substituting these into the PB Eq. (|1.8[) we find: 

2 / o'T' \ 1 / T n \ 1 s 



2 V T 



m + n 



= 0. (3.2) 



The far field limit z — > oo corresponds to s = 0. Around 
this point, T(s) can be expanded into a Taylor series: 



T( S ) 



OO 

fe=0 



1 + C\ s ■ 



c 2 s 2 + 



(3.3) 



The zeroth order coefficient cq = 1 because as z — > oo, 
s -» 0, * ->• 0, and T = e* -» 1. Substituting the expan- 
sion Eq. (|3.3|) into Eq. (|3.2j) and compare the coefficients 
of s fc order by order, we easily obtain relations between 
Cfc's, which allow us to solve all c^'s in terms of c\. 



(3.4) 



The numbers Ck can be calculated straightforwardly us- 
ing Wolfram Mathcmatica up to arbitrary order. In 
this work we only keep twelve terms for (m : — l)(m = 
1,2,3,4) electrolytes, and sixteen terms for (3 : —2) and 
(4 : —3) cases. Some of these coefficients Ck are tabu- 
lated in Table HU The only remaining coefficient c\ will 
be determined by matching the far field expansion with 
the near field expansion, which, by construction, is due 
to an infinitely charged plate at the origin. 

Combining Eq. (|3.3|) and Eq. (|3.1[) we can also obtain 
the far field expansion for the potential To the leading 
order we have: 



#(z) = c ie - z +0(e~ 2z ). 



(3.5) 



Comparing this with Eq. (| 1 . 5|) and Eq. (|1.11[) (also notic- 
ing the relation between </> and VP) we see that 



ci = rj R (±oo), 



(3.6) 



i.e., ci is precisely the dimensionless renormalized dimen- 
sionless surface charge density for an infinitely charged 
plate. 



To employ asymptotic matching of the far field expan- 
sions and the near field one, we will have to evaluate 
the truncated far field expansion at z* ~ 1. However, It 
turns out that, for many cases of (m,— n), the far field 
expansion does not converge near z* ~ 1. [l9[ We use 
the well known Shanks transformation (TBI to improve 
convergence for the far field expansion of T(s). In Ap- 
pendix we give a brief introduction to this method. 
The useful Shanks transformation is of second order. 
Consider a power series A = E„ a n x n . The partial sums 
is A n = Efe<« a nX n , its second order Shank transform is 
given by 



S 2 (A n ) 



A n -2 


A n -1 


A n 


AA n ^ 2 




AA n 




AA n 


AAn+1 


1 


1 


1 


AA n ^ 2 


AA n -l 


AA n 


A4-1 


AA n 


AAn+1 



(3.7) 



where AA n = A n — A n -± = a n x n . Replacing of A n by 
S2(A n ) may lead to substantial improvement of conver- 
gence. 

Implementation of the second order Shanks transform 
is straightforward using Wolfram Mathematica. For the 
cases of 1 : — 1 and 2 : — 1, Shanks transform yields very 
simple result: 



- Shan 



k(*) 



(4+cisf 
(4-ci s)' 2 ' 



(6~c lS ) 2 



1 : -1; 

2 : -1. 



(3.8) 
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(~2 


C3 




C5 


C6 


Or 


C8 


<"<> 


CIO 


Cll 


C12 


1 


-1 


1 


3 


1 


r > 


3 


7 


1 


9 


5 


11 


3 


2 


16 


16 


25(1 


512 


4096 


2048 


65536 


131072 


1048576 


1048576 


2 


-1 


1 


1 


1 


5 


1 


7 


1 


1 


5 


11 


1 


3 


12 


54 


1296 


1296 


40656 


34992 


186624 


5038848 


60466176 


30233088 


3 


-1 


1 


1 


5 


269 


125 


21469 


25997 


76699 


10511275 


216900209 


693028429 


6 


16 


432 


20736 


13824 


2985984 


4478976 


15925248 


61917364224 


61917364224 


227030335488 


4 


-1 





1 

8 


1 

10 


33 
320 


11 
100 


313 
2560 


1567 
11200 


83853 
512000 


1749 
8960 


33847223 
143360000 


22792479 
78848000 


3 


-2 


1 


1 


7 


49 


53 


941 


1 


673 


4549 


109487 


252281 


3 


6 


135 


1620 


8100 


145800 


54675 


364500 


7873200 


147622500 


590490000 


4 


-3 


1 


7 


11 


605 


935 


146399 


17149 


4797881 


6268525 


2229207493 


703428293 


3 


24 


108 


5184 


36288 


2612736 


13716864 


146313216 


564350976 


94810963968 


47405481984 



TABLE II: Some coefficients Ck of far field expansions. 



which exactly solve Eq. (|3.2[) . They correspond to exact 
solutions to the original nonlinear PB equation Eq. (jl.8l) : 



*(*) 



2in g +cie :;i 

(4 — cie z ) ' 



, 36+24cie~"+cje~ 
111 (6-cie-*) 2 



(3.9) 



The parameter c\ can be determined by matching with 
the boundary conditions on the plate. It is interesting 
to note that in the far field expansions c\ always appear 
together with e~ z , hence change of c\ amounts to trans- 
lation of the coordinate z. 

For other valences m : —n, the Shanks transforms of 
the far field expansions do not yield exact solutions. Nev- 
ertheless, they remarkably improve the quality of approx- 
imation. We substitute the far field expansions (with c\ 
set to unity) as well as their Shanks transformed versions 
into the nonlinear PB equation Eq. (|3.2[) . Since they are 
not exact solutions, the right hand sides do not identi- 
cally vanish, but measure the quality of the truncated far 
field expansions. As shown in Fig. [TJ the Shanks trans- 
form reduces the error of the far field expansion by about 
three orders of magnitude for the case of (3 : —1). Similar 
improvements are also achieved for other cases. 



IV. ASYMPTOTIC MATCHING 



We now require that the far field expansions equal to 
the near field expansions (both of them are, of course, 
approximate) at some intermediate point z*: 



*fa P r P (^Vl) 



* app (z*) 



Note that only the far field expansion depends on c\. 
Solving this equation determines the only parameter in 
the far field expansion, i.e. the rcnormalizcd charge den- 
sity c\ = ?7i{(±oo). The matching condition Eq. (|4.1[) 
is numerically solved for c\ using Mathematica. The 




FIG. I: The error of truncated far field expansion as an ap- 
proximate solution to the nonlinear PB equation Eq. (|3.2[) , as 
well as that of the Shanks transformed version, for the case 
of (3 : — 1). The Shanks transform reduces the error by about 
three orders of magnitude. 




FIG. 2: Matching of near field and far field expansions for a 
negatively charged plate in 3 : —1 electrolyte. Top: far-field 
expansion, Bottom: near field expansion. Evidently, the far 
field expansion is excellent even around tenth of Debye length. 



(4.1) solution however can only be an approximation to the 



exact value cf x , since both near field and far field expan- 
sions are approximate. The errors of these expansions 
can be estimated by the difference between the current 
order approximation and the next order approximation. 
Via some simple analysis, to be detailed in Appendix IB! 
we derive the following approximate bound for the error 
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electrolyte 


plate 


z* 


Afar 


ticar 


(12 


ci = ?7ii(±oo) 


5ci 


1 : -1 


+ 


I 

1 


12 


/ 

1 


/ 

1 


4 


exact 




1 


12 


1 


1 


-4 


exact 


m : —m 


+ 


1 

1 


12 


1 

1 


1 

1 


4/m 


exact 






1 


12 


1 


1 


— 4/m 


exact 


2 : -1 


+ 


I 

1 


12 


1 

1 


I 

1 


6 


exact 




1 


12 


1 


1 


-6(2 -V3) 


exact 


1 : -2 


+ 


I 

1 


12 


I 

I 


1 

1 


6(2 - \/3) 


exact 






1 


12 


1 


1 


-6 


exact 


2m : — rn 


+ 


I 

1 


12 


I 

1 


I 

1 


6/ttj. 


exact 






1 


12 


1 


1 


-6(2 - \/3)/m 


exact 


3 : -1 


+ 


2.0 


12 


li 


1/18 


8.707001 


8 x 10 -6 






0.5 


12 


10 


-1/6 


-0.9938796 


2 x 10~ 7 


1 : -3 


+ 


0.5 


12 


10 


1/6 


0.9938796 


2 x 10" 7 






2.0 


12 


11 


-1/18 


-8.707001 


8 x 10~ 6 


4 : -1 


+ 


2.5 


12 


11 


1/24 

/ 


12.3142 


0.0002 






0.4 


12 


10 


-1/6 


-0.717417 


2 x 10 -6 


3 : -2 


+ 


0.7 


16 


13 


1/18 


2.45953455 


1 x 10~ 8 






0.6 


16 


9 


-1/12 


-1.1542225 


3 x 10~ 7 


4 : -3 


+ 


1.0 


16 


13 


1/24 


1.53280486 


1 x 10~ 8 






0.6 


16 


7 


-1/18 


-0.89788695 


7 x 10~ 8 



TABLE III: Saturated values of renormalized surface charge density tir(±oo) for different electrolytes. iVf a r and A^noar are 
numbers of terms kept in the far field and near field expansions. Two identities that can be used to generalize the above 
results to other cases: r/^'~ n (oo) = — ?^' _m (— oo). and ?7^'~ n (oo) = p _1 ^' _n (oo), where m = pm,n = pn. 



Sci = Ci 



\6a\ < 



\S^(z*)\ + \S^(z*, Cl )\ 



(4.2) 



dci 



We systematically vary the matching point z* to mini- 
mize this error bound. The optimal z* , the calculated 
coefficient c\, and the corresponding error bound for all 
cases are shown on the Table. Mil Except for the case of 
positively charged plate in 4 : — 1 electrolyte, the errors 
are all with in orders of 10~ 6 . Also shown there are the 
number of terms used in the near field and far field ex- 
pansions. In Fig. [5] we show both near field and far field 
expansions for the case of a negatively charged plate in 
3 : — 1 electrolyte. It can be seen there that the Shank 
transformed far field expansion is excellent even if z* is 
within a Debye length to the plate. 

Note that our asymptotic matching is different from 
the one discussed in standard textbooks of asymptotic 
analysis. There, it is required that there is a control- 
ling small parameter e, and a whole intermediate region 
where both far field and near field expansions become 
exact as e — > 0. It is then guaranteed that in this limit, 



the matching becomes exact. In our problem, there is 
no small parameter whatsoever. The precision of both 
far field and near field expansions are guaranteed by the 
high order of expansions, augmented by Shanks transfor- 
mation. Evidently, analysis with these high order expan- 
sions would be extremely time consuming without mod- 
ern computer softwares for symbolic calculation, such as 
Wolfram Mathematica or Maple. 



RENORMALIZED CHARGE DENSITY OF A 
HIGHLY CHARGED PLATE 



Using the calculated coefficient c\ = 7yij(±oo), we can 
now completely determine the far field asymptotic expan- 
sion of the potential ^(z). Recall the near field expansion 
of *&(z) exhibits a logarithmic singularity at z = 0, see 
Eqs. (|2.2I2.4[) . which corresponds to an infinitely charged 
plate located at the origin. What we really want is how- 
ever the asymptotics of a finitely charged plate, say, with 
dimensionless surface charge density r\. We can always 
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Electrolytes Plate Asymptotic expansion of r)a(jf) 



3 : -2 



7=22 (exact) 



, fi _12 , 12_ 4 _ 4 , 4 , 4 _ 4 , 4 , 92 _ 124 , n / 

- 6(V3-2) + ^ + £(^ + ^ + ^ + ^ + ^ + o ( ^ 

, o 707Q1 (l — 1 + 2 — 0-888888 _ 0.222225 . 0.444446 _ 0.0493818 _ 0.197532 \ , q ('•j-SN 

- -0.993871 (l - ^ + ^ + _ <L™ _ M~ + OM44T04 _ a™) + Q ^,-ir/a) 

, ^ 3141 (l_2± 2 _ 0-999991 _ 0. 0000187721 , 0.300013 _ 0.0999981 _ 0.0785797 \ , q („-8\ 

' \ V V 7 ?P V® i 77 / V ' ' 



717419 (1 — + 1 + ".0625227 _ 0.0401C72 _ 0.0390738 , 0.0200836 _ 0.0148545 

lil sFTP FTP hi 7 / 2 FTP 3 M 9 / 2 FTP 5 



+ o(V 11/2 ) 



, 2 46097 ( 1 — - + 2 — °- u55 J 251 _ 0-0695749 i 0.0250130 _ 0.00634138 _ 0.00260921 \ . q (^-S 1 ) 

— —1 154222 ( 1 - h 2 4- o o 2 5 21 92 _ 0.0415042 _ 0.0133842 . 0.000454105 . 0.00892278 \ , q (\ rl \- ( >) 

\ 3|t7 ' 9 1 1; 1 2 FTP FTP 3 |r;| 13 / 3 FTP 5 |r;| 16 / 3 J V ' / 



0.0123480 0.0164594 , 0.00342930 0.00085279 , 0.00144768 \ , r> ^~ 2 0/3\ 

-r ^'"^j^t^ ? ^ + ? ^773— + ^ )+0[r) ) 

— 897888 ( 1 - h 1 4- 0-00694196 _ 0.0112835 _ 0.00255865 . 0.000434214 , 0.00127932 \ . q n i 



"The exact result is given implicitly by ri = 3S vrXvr+6) 

6 (6- vr )( Vr + 24. vr +36) 



TABLE IV: Asymptotic expansions of the renormalized dimensionless surface charge densities 77^(77) in terms of the bare 
dimensionless density 77, for different electrolytes. Using Eq. p. 151 [TTTS)) . we can generalize these to other cases as well. 




FIG. 3: Renormalized surface charge densities as functions of bare surface charge densities. Top: positively charged plate. 
Bottom: negatively charged plate. 



set up the coordinate system such that the 77 plate is lo- 
cated at zq (cf Eq. (|1.9|l ). such that if we extrapolate the 
solution &(z) to the region z < zq, there is a singularity 
at z = 0. This means that the potential $(z) is identical 
to the one have analyzed in previously. The location Zq 
of the plate now becomes a function of 77, through the 
Neumann boundary condition Eq. (|1.9[) . Using the near 
field expansions for ty, the boundary condition Eq. (|1.9|) 
now becomes 

- "^near^u) = V- (5-1) 



The left hand side can be expanded into an asymptotic 
series in z 0l using the near field expansion discussed in 
Sec. HH The series can be inverted to yield z as an asymp- 
totic expansion in terms of 77 = frjg/^, which is also 
small for strongly charged surface. The results are how- 
ever to complicated to be shown explicitly. 

On the other hand, the leading order far field asymp- 
totics can be written as 

~ ^(±00) e~ z = Vr (±oo) e~ z ° e~ Sz , (5.2) 

where 8z = z — zq is the distance from the plate to the 
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field point. Therefore by the definition Eq. (jl.5p . the 
renormalized dimensionless surface charge density of the 
plate with bare surface charge 77 is given by 

m(v)=VR(±°o)e- z °. (5.3) 

Now expanding e~ z ° into asymptotic series of 77, we 
finally obtain the renormalized surface charge density 
Vniv) as an asymptotic series in terms of the bare sur- 
face charge density 77. The results are shown in Table 
IIVI and Fig. [3] Even though these results are derived for 
strongly charged surface, they are quantitatively accu- 
rate even for weakly charged surfaces. For example, for 
the case of 7/ = £db/h = 5, the relative error for 1 : — 1 
electrolyte is of order of 10~ 6 . 



VI. CONCLUDING REMARKS 

In this work, we have developed systematic near field 
and far field expansions for the electrostatic potential 
generated by a strongly charged plate inside a generic 
charge-asymmetric electrolyte. Using these expansions, 
we have derived a perturbative expansion for the renor- 
malized surface charge density in terms of the bare sur- 
face charge density. Our methods can be immediately 
generalized to charged plates immersed in mixed elec- 
trolytes. 

The authors thanks Shanghai Jiao Tong University and 
NSFC (grant numbers 11174196, 91130012) for financial 
support. 
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Appendix A: Basics of Shanks Transform 

In this appendix, we present a heuristic discussion of 
Shanks transform. The readers are referred to the clas- 
sical monograph by Bender and Orszag [l5[ for details. 
Originally Shanks transformation was invented to speed 
up convergence of slowly converging series. Its most 
powerful application is however to sum up divergent se- 
ries and yield finite result. For simple geometric series, 
Shanks transformation is equivalent to the traditional 
procedure, where one first sum up the series within its 
radius of convergence, and then analytically continuate 
beyond this domain. 

Consider a simple geometric series 




which is convergent for \z\ < z\. Within this domain 
of convergence, the series can be summed up to yield a 
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simple rational function with one simple pole at z = Z\ 

a 



A 



l-z/zi 

Let us define the remainder R n as 
R n = A A n . 



(A2) 



(A3) 



which is the error if we use the partial sum A n to approx- 
imate the original series. For the simple geometric series 
Eq. (|A1|) the remainder has the following simple form 



Rri 



1 — z/zi \zi 



(A4) 



which is also geometric in n. The remainder goes to zero 
as long as the original series converges. 

Suppose we are interested in the domain \z\ > 1, i.e. 
outside the radius of convergence of A, and suppose that 
for some reason, we are only able to calculate the partial 
sums A n = a z k only up to some finite order n. The 
conventional approach fails here because we can not sum 
up the series to infinite order for small z. On the other 
hand, the partial sum A n is not a good approximation 
to the exact because the error (i.e. the remainder R n ) 
diverges in the domain of interest. So the question is 
how we can obtain a good approximation to the series 
outside its domain of convergence, with the knowledge of 
only finite terms of the series? 

For the geometric series Eq. (|A1|) . we have the follow- 
ing relations for arbitrary integer n: 



A n +i 


= A — 


az n+1 , 


(A5) 


A n 


= A - 


a z n , 


(A6) 


A n -i 


= A - 


az n ~\ 


(A7) 



Solving these three equations, we can express A in 
terms of three parameters A n+ i,A n ,A n ^\ (details ig- 
nored here): 



.4 




(A8) 



where AA n = A n+ \ — A n . Therefore knowing that the 
series is geometric, we only need three terms of partial 
sums to obtain the exact sum of the series. 

The Shanks transformation is a nonlinear transforma- 
tion acting on the partial sums Si(A n ) (The subscript 1 
means it is the first order Shanks transformation). 



Si (AO 



A n +iA r , 



A'i. 



A 



71+1 



An-l — 2A n 



(A9) 



Note however, in computing Si(A n ), we actually need 
three partial sums A n+ i, A n , A n -\. Strictly speaking, 



therefore, Shanks transformation acts on the whole se- 
quence of partial sums, and return a new sequence of 
partial sums. It is defined such that it returns the exact 
A if the original series is geometric. 

The utility of Shanks transformation can be best un- 
derstood by looking at its effect on the remainder. Let 
A n = A — R n in the above equation, where R n is the 
remainder of the original series. Further let R^ be the 
remainder of the Shanks transformed scries: 



A - Rt 



Si(A n 

We easily see the following 

Rn+lRn-l 



R„ = 



R7, 



Rn+l + Rn-1 — 2R n 



(A10) 



(All) 



If the original series is geometric, then R n has the form 
of Eq. (|A4[) . and R^ simply vanishes, indicating that the 
transformed partial sum is actually the exact result. 

For a more general series, let us assume that we can 
separate the dominant part of the remainder which is 
again gemoetric: 



Rn = R n + a z r ' 



(A12) 



where R n /z n — > as n — > 00. Our analysis below applies 
regardless of the convergence of the remainders R n ,R n . 
Using Eq. (|All|) . the Shanks transformed remainder then 
becomes: 



R', 



a(z 
zR n 



a(zR n -i + z 1 R n+ i - 2R n ) 

- z- 1 - 2) + (R n+1 + it,,-! - 2R n )/z n 
-1 + z~ 1 R n +\ — 2R n 



(z + z- 1 ~2) 



(A13) 



In the limit n — > 00, therefore, the dominant part z n no 
longer appear in the transformed version of the remain- 
der. As a consequence the transformed remainder be- 
comes much smaller than the original remainder. In par- 
ticular, if z > 1 but R n — > 0, the original series diverges 
Rn — > 00 but the Shanks transformed series converges, 
i.e. R% — > 0. This shows how Shanks transformation can 
be used to achieve the purpose of analytic continuation 
with only finite number of terms of the series. Note that 
for the case we studied here, the summation Eq. (|A2[) 
contains a simple pole at z = 1. The effect of Shanks 
transformation can also be understood as removing this 
singularity from the original series. It is easy to see that 
the Shanks transformation recovers the exacts from Tay- 
lor series of all rational functions with only one pole, in 
the form of P(z)/(z — zq), where P(z) is a polynomial of 
arbitrary (but finite) order. 

The first order Shanks transformation discussed above 
does not work if the remainders contain second order 
poles. To deal with these series, we have to introduce 
second order Shanks transformation, which removes two 
poles at once. In this case, it can be shown that the exact 
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sum A can be calculated using five terms of the partial 
sums: A n+2 ,A n+ i,A n ,A n _i,A n _ 2 - 



A 





An-i 




2 


AA n 




AA n 


-2 


AA n +i 


AA n 


AA n 


-1 


1 


1 


1 




AA n 




AA n 


-2 


A4+i 


AA n 


AA n 


-1 



(A14) 



For a general series, the second order Shanks transfor- 
mation is then defined as 



S 2 (A r , 



A n 


A n -1 


A n - 


2 


AA n 


AA n -! 


AA n 


-2 


AA n+ i 


AA n 


AA n 


-1 


1 


1 


1 




AA n 


AAn-1 


AA n 


-2 


AA n+ i 


AA n 


AA n 


-1 



(A15) 



This second order transformation returns the exacts 
when applying to the Taylor expansions of rational func- 
tions with two poles. 



Appendix B: Estimation of Error 5ci 

To estimate the error of Ci , we note that if we use the 
exact near field and far field expansions in Eq. (|4.1[) . we 
would have found the exact value for the coefficient cf x , 
independent of the matching point z* . Mathematically 
we have: 



where the superscript ex stands for "exact" . By contrast, 
because of the approximate nature of the expansions we 
use, the coefficient c\ we find is different from its exact 
value. Let c\ = c° x -I- 8c\ 1 where 8c\ is the error, we can 
rewrite Eq. (|4. 1 [) in the following form: 



*l™(z*,c? + 5c) = *%l(z*). (B2) 

Expanding the left hand side to the linear order of Sc\, 
we have: 

(B3) 

The derivative of \&f PP with respect to c\ can be easily 
calculated using Mathematica. We further introduce the 
error of the far field and near field expansions (at the 
matching point z*) as following: 

6* iaI (z*,c) = ^ a p r p (z*, Cl )-* far (2*, Cl ), (B4a) 

**ne«CO = - *near(«*). (B4b) 

These errors can be estimated by the difference between 
the approximation at the current order and that at the 
next order, which can be conveniently calculated using 
Mathematica. Using Eqs. (|B1[ IB31 IB4[) . we can estimate 
the upper bound of the error 5c as 



\Sc\ < 



|**£ p (** )Cl )| 



(B5) 



dci 



Appendix C: Coefficients of Near Field Expansions 



*£(«*,cf) = *S ar (z*) ) 



(Bl) 
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